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Abstract 

A model based on a nonrelativistic potential supplemented by spin 
dependent terms of the Breit-Fermi type is described. This model re- 
produces the known masses of the bb quarkonia within lMeV and the 
known masses of the cc quarkonia within 5MeV. Many predictions are 
made, some of them, e.g. for the fine splittings of the high L states of 
the bb quarkonia and of the mass spectrum of the be quarkonia, sig- 
nificantly different from previous predictions made by other authors. 



1 Introduction 

For the hydrogen atom the nonrelativistic Schrodinger equation yields a good 
first approximation for the energy levels. The spin effects are responsible only 
for the fine and hyperfine structure of the levels. In heavy quarkonia, as well, 
the spin effects yield the fine and hyperfine splittings, but there are two im- 
portant complications. In the hydrogen atom fine splittings are negligible 
compared to the energy differences between the nonrelativistic levels and 
hyperfine splittings are negligible compared to the fine splittings. In heavy 
quarkonia the fine and hyperfine splittings are of the same order. They are 
smaller than the energy intervals between the unsplit levels, but certainly 
not negligible. Therefore, the comparison of nonrelativistic calculations with 
experiment is ambiguous. In the hydrogen atom the velocity of the electron 
in the ground state is about 0.01c, thus the nonrelativistic approximation is a 
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well justified starting point. In heavy quarkonia the quarks are much faster. 
The root mean square velocity is, in the ground state of the bb system, about 
0.3c and in the ground state of the cc system about 0.5c. The c-quark in the 
be system is probably even faster than that in the cc system. It may seem, 
therefore, that the nonrelativistic approximation is suitable only for crude 
estimates. Nevertheless, many authors have found that nonrelativistic mod- 
els give amazingly successful predictions for the spin averaged energy levels 
of heavy quarkonia. Including perturbatively spin effects one can also repro- 
duce very satisfactorily the fine and hyperfine structure. Many references to 
work of this kind can be found in the review paper J^j. A simple-minded 
interpretation is that most of the relativistic effects is somehow taken into 
account by fixing the parameters of a nonrelativistic model to fit the exper- 
imental data - they just renormalize the parameters of the nonrelativistic 
potential. It could be argued that instead of the poorly understood po- 
tential models it is better to use QCD sum rules (cf. e.g. 0) or lattice 
methods (cf e.g. [Q). At present, however, these methods have rather large 
uncertainties. In the present paper we describe a model 0,0, which gives 
a remarkably good description of the mass spectra below the strong decay 
thresholds, i.e. below 10.558 GeV for the bb quarkonia and below 3.729GeV 
for the cc quarkonia. This makes it plausible that also its predictions for the 
be quarkonia below 7.144 GeV (7.189 GeV) will not be far off the mark. The 
mass limit in brackets is applicable to quarkonia, which, because of angu- 
lar momentum and parity conservation, cannot decay into two pseudoscalar 
mesons. In these ranges there are 34 mass states for the bb quarkonia (out 
of that nine are known), 16 mass states for the be quarkonia and 8 states for 
the cc quarkonia (out of that six are known) . We consider as known only the 
mass states included as firmly established in the 1986 Particle Data Group 
Tables §. 

2 Nonrelativistic potential 

Our model contains six free parameters - the quark masses 

m b = 4.8030GeV, m c = 1.3959GeV, (1) 
the coupling constant a s at the m c mass scale 
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a s (m c ) = 0.3376 
and the three parameters of the nonrelativistic potential 



(2) 



a = 0.32525, b = 0.70638GeV5, c = 0.78891GeV. (3) 

Of course we do not claim to have fixed the quark masses with an accuracy of 
four digits etc. All these digits are necessary, however, in order to reproduce 
our results. The potential is spherically symmetric and reads 

a 

V (r) = mQ + m-Q — c + b\fr . (4) 

The eigenvalues of the corresponding Schrodinger equation are interpreted 
as the centres of mass of the spin one multiplets. For the be quarkonia this 
corresponds to the case, when the singlet-triplet mixing is switched off. 

Comparison with experiment for the bb quarkonia can be made in five 
cases. The agreement is within 1 MeV. A rough estimate of the relativistic 
correction gives about 60 MeV. Thus, almost all this correction must reduce 
to a redefinition of the parameters of the potential. For the cc quarkonia 
comparison with experiment can be made in three cases. The dicrepancies 
are within 4 MeV. Since the relativistic corrections here should be much 
larger than in the previous case, the good agreement is even more surprising 
than before. It is instructive to compare our results with those of Gupta and 
Johnson |7|]. They are using a very different model, but their fit is about 
as good as ours. In fact it is somewhat better for the cc system, though 
somewhat worse for the bb system. In spite of the agreement within a few 
MeV between the two models for the bb and cc systems, the predictions for 
the be quarkonia differ significantly. We expect higher masses, by about 40 
MeV for the S states and by 16 MeV for the P multiplet. Thus, a good fit for 
the bb and cc is compatible with a poor fit for the be quarkonia. Of course, 
we do not know yet, whose fit is poor, but both fits cannot be good. 

3 Hyperfine splittings 

The hyperfine splittings are usually described using the Breit-Fermi interac- 
tion, which in the present case takes the form (cf. e.g. f|) 
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In the first approximation of the Rayleigh-Schrodinger perturbation theory 
the corresponding energy splitting is 

In this approximation, the splitting occurs only for the 5 states, because 
for the states with L ^ the wave function ip vanishes, when the distance 
between the two quarks goes to zero. The measured splitting of the 15* state 
of the cc system yields the value of a s (2). Once this is chosen, it is possible 
to predict any other hyperfine splitting, but there is no data to compare 
the predictions with. Some indirect support is obtained by calculating the 
coupling constant a s at the mass of the Z° — the result a s (M z o) = 0.115 
is very acceptable — and from the calculation of the leptonic decay widths 
of heavy quarkonia |J. It is instructive to compare our prediction for the 
hyperfine splitting of the IS state in the bb system (56.7 MeV) with the 
predictions of the sum rules @ — 63 ±30 MeV — and of the lattice approach 
[||] - 60 MeV with a large, but unspecified, error. Our agreement with 
the mean values is good, but in view of the large uncertainties of the more 
rigorous approaches, it is premature to draw any conclusions. 



4 Fine splittings and mixing 

The interaction terms responsible for the fine splittings are (cf. e.g. ||, |T0|| ) 
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In these formulae the only undefined quantity is the coupling constant a s , 
which should be running, i.e. depending on r. We propose the following 
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procedure to define this function. The general formula is modelled on the 
well-known one-loop expression for a s (p) with the replacement of the momen- 
tum scale p by K We introduce tildes above a s and A, when considered as 
functions of r, in order to avoid confusion with the corresponding quantities 
considered as functions of momentum. Thus 



12- (A<"'>)V (iq) 
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Here the number of flavours rif equals three for r > m" 1 , equals four for 
< r < m" 1 and equals five for r < . The region r > m" 1 is of little 
interest for the present problem, therefore, we keep rif = 3 also for large 
values of r. The numerator is introduced in order to compensate the zero of 
the denominator at r = . Its exact shape is of little importance. The 

parameter A^ 4 ^ is obtained from the two conditions 



in, 



a s (m c ), (11) 



a s (—) = a s (m b ). (12) 
\m b J 

These two equations give slightly different values of the parameter A( 4 ) and 
we take the geometrical mean of the two results. However, the two solutions 
are so close to each other, that taking the arithmetic mean instead of the 
geometric one would make no difference at our level of precision. Then, A (3) 
and A^ 5 ) are calculated from the conditions that the function a s (r) should be 
continuous at r = m" 1 and at r = m^ 1 . Interpreting a s in the formulae for 
the potentials Va, Vls an d Vt as the function a s (r), we obtain the explicit 
form of the operator responsible for the fine splittings of the mass levels of 
the heavy quarkonia. 

When calculating the splittings of the bb quarkonia, we find very good 
agreement with experiment for the splittings of the IP and 2P states. The 
accuracies of our results are all within lMeV. Most other models, however, are 
also successful in these predictions, typical errors being within 5 MeV. Also 
for the 3P states, where there is no data to compare with, our predictions 
are close to those of other authors. A qualitative difference between the 
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predictions appears, however, for the higher L states. For instance, Kwong 
and Rosner |TTJ predict a splitting of about 1 MeV for the F-states and 
almost no splitting for the G-states, while we find a splitting of about 13 
MeV for the IF state and of about 10 MeV for the 1G state. Thus, the 
high L states will be of great interest for the comparison of models. For 
the cc systems our predictions for the splitting of the IP state agree with 
experiment within about 5 MeV. For other states we are in rough agreement 
with other models, except for the 3 Di state, where we predict a down shift 
of about 20 MeV, while other models find much smaller shifts. 

For the be system one can calculate the singlet-triplet mixing due to the 
operator Va- These angles are very sensitive to the details of the models 
and comparing the results of our paper and of the papers || and |IJ one 



sees that no two sets of predictions agree. According to our calculation the 
sines of the mixing angles for the IP, 2P and ID states are respectively: 
0.374, 0.385 and 0.244. Eichten and Quigg § find almost no mixing; for the 



IP state, while Gershtein and collaborators [1C] find the biggest mixing for 
the ID state. 

Another interesting mixing effect is due to the operator Vp- This involves 
pairs of states differing by two units in orbital angular momentum. The 
mixing angles are very small — less than 10~ 3 — and the effect on the mass 
values is negligible. This effect, however, enhances considerably the leptonic 
decays of the quarkonia with angular momenta L > 2. 
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